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Abstract. We investigate the effect of geometrical frustration on the H-T phase 
diagrams of the classical Heisenberg antiferromagnets on triangular and kagome 
lattices. The phase diagrams for the two models are obtained from large scale 
Monte Carlo simulations. For the kagome antiferromagnet thermal fluctuations 
are unable to lift completely degeneracy and stabilize translationally disordered 
multipolar phases. We find a substantial difference in the temperature scales of 
the order by disorder effect related to different degeneracy of the low- and the 
high-field classical ground states in the kagome antiferromagnet. In the low-field 
regime, the Kosterlitz-Thouless transition into a spin-nematic phase is produced 
by unbinding of half-quantum vortices. 
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1. Introduction 

One of the important issues in the field of highly frustrated magnets is the universality 
of the order by disorder phenomenon [TJ [2]. If the 'accidental' degeneracy of 
a frustrated spin system is described by a few continuous degrees of freedom 
[H [3j |H O [7] or by an infinite but non-extensive number of such parameters [9] 
the order by disorder mechanism is known to stabilize a magnetically ordered state 
at low temperatures. In the case of a macroscopic number of zero-energy modes the 
outcome is less universal and only a partial lifting of the degeneracy, if any at all, may 
take place QUI [TTJ [H] . 

In this work we study numerically the entropic order by disorder selection 
produced by thermal fluctuations in two classical antiferromagnets on a triangular and 
a kagome lattice in external magnetic field. The external field changes continuously 
the classical ground-states manifold of a frustrated magnet providing an experimental 
tool to control and to modify the effect of fluctuations [T31 [2] ■ The magnetic phase 
diagrams for the two models have been discussed qualitatively in the earlier works 
p~5l [To] . Here we present results of large-scale Monte Carlo simulations, which 
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Figure 1. Triangular (a) and kagome (b) lattices. Equilibrium three-sublatticc 
configurations in different field regimes: (c) for < H < |H S , (d) for H = \H 3 , 
(c) for \H S < H < H s and (f) for H > H s . 

allow precise mapping of the transition boundaries and clarify nature of the low- 
temperature phases. In addition to providing new insights into the order by disorder 
effect in classical frustrated magnets, our results open a possibility of quantitative 
comparison with the experimental phase diagrams of recently discovered triangular 
antiferromagnets RbFe(Mo0 4 ) 2 Q2H3U and Rb 4 Mn(Mo0 4 ) 3 [31]. 

The nearest-neighbor exchange Hamiltonian with a Zeeman term is given by 

(ij) 

where Si are three-component classical unit vectors. The model (pj with an 
antiferromagnetic exchange J > is studied for a triangular and a kagome lattice, see 
figures [31 a) and (b). The common building block of the two lattices, a triangular 
plaquette, determines frustration at the shortest length-scale. In zero field the 
antiferromagnetic bonds on a single spin-triangle cannot be simultaneously satisfied 
and magnetic moments form a 120° spin-structure. For the triangular antiferromagnet 
this leads to an ordered magnetic state with the ordering wavevector Q = (47r/3, 0) 
and a trivial six-fold degeneracy related to permutations of the spin triad. In contrast, 
the kagome-lattice model remains infinitely degenerate. The coplanar ground-states 
of the kagome antiferromagnet can be mapped onto different colouring patterns of 
an exactly solvable three-colour problem [20l [21]. The number of colourings grows 
exponentially with the number of sites N as 

W « 1.1347^ . (2) 

In a finite magnetic field the classical ground states may be derived by rewriting 
the energy flj) as a sum over triangles [131 HI]- The classical energy reaches the 
minimal value for the following magnetization of every plaquette: 

E s * = 3 f-^ ( 3 ) 

where H s stands for the saturation field with H s = 9 J and 6 J for the triangular and 
the kagome antiferromagnet, respectively. The constraint ([3]) leaves three parameters, 
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instead of one, for the choice of the spin triad Si at H < H s . Physically, extra 
degrees of freedom correspond to an isotropic susceptibility of a spin triangle. Thus, 
an external magnetic field produces new degeneracy in the case of the triangular 
antiferromagnet and enhances the existing degeneracy for the kagome model. 

The entropic degeneracy lifting in the classical triangular antiferromagnet was 
first studied by Kawamura, who calculated the free-energy contribution of the 
harmonic excitations [3]. Thermal fluctuations lift the degeneracy and stabilize in 
the interval < H < H s one collinear and two coplanar states. These are shown 
in figures (Hc)-(e) and will be referred to in the following as the coplanar Y-state, 
< H < |i? S) the coplanar V-state, |if s < H < H s , and the collinear up-up-down 
(uud) state, H ~ j$H s . 

The harmonic fluctuations in the kagome antiferromagnet are dominated by zero- 
energy modes, which contribute most to the free-energy decrease at low temperatures 
[2D] . Their number for different classical ground states (J3{ can be counted by using 
purely geometrical arguments. As a result, one ends up with the same sequence of 
three-sublattice configurations in magnetic field, figures HJc)-(e) [16]. In addition, the 
zero-energy modes favour the \/3x \/3 periodic pattern for the V-state, whereas no such 
selection takes place for the Y- and for the irad-state. Note, also that two sublattices 
become identical for uud and V configurations. Consequently, the total number of 
these states is smaller than ([2J and is given by the number of dimer coverings of the 
dual hexagonal lattice [23] [16] : 

1.1137^. (4) 

Lower degeneracy in the high-field region \H S < H < H s may significantly affect the 
fate of the order by disorder selection. 

In the following we present the Monte Carlo results for the phase diagrams of the 
triangular (section [2]) and of the kagome antiferromagnet (section l3.2p . In section I3TT1 
we also give a brief symmetry analysis of various multipolar states in zero and finite 
magnetic fields. The obtained results are summarized in section 2] 



2. Triangular antiferromagnet 

Monte Carlo simulations have been performed using the standard Metropolis algorithm 
in combination with the microcanonical over-relaxation steps, see [H] for further 
details. Periodic boundary conditions were implemented for N = L x L site clusters 
with the linear size L up to 192. At every temperature/magnetic field we discarded 
5 x 10 4 Monte Carlo steps (MCS) for initial relaxation and data were collected during 
subsequent 10 6 MCS. The error bars were estimated from 10-20 independent runs 
initialized by different random numbers. Instantaneous values of the antiferromagnetic 
order parameter 

m Q = ^E s * e ~ iQri ' Q = (!^o) (5) 

i 

have been used to measure longitudinal and transverse components of the staggered 
spin susceptibility and the corresponding Binder cumulants: 

N „ ((™?>) 4 ) 

XQ=y<(™q) 2 >, ^Q = ^^. (6) 

where a — z, _L and iUq — (iBq + TOq ) 1 / 2 . 
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Figure 2. Transition between the paramagnetic phase and the uud-state in the 
triangular antiferromagnet at H/J = 1: the Binder cumulant Uq (upper panel) 
and the staggered susceptibility Xq (lower panel) versus temperature for different 
lattice sizes L. 



Following the standard procedure, the second-order transition may be located 
by the crossing point of the Binder cumulant Ul(T) measured for different clusters. 
We illustrate this method in the top panel of figure [2] for the transition between the 
paramagnetic state and the uud state at H/J = 1. The alternative approach is to use 
the order parameter susceptibility provided the critical exponent r\ is known precisely. 
In the critical region the susceptibility scales as 

x = L 2 -Vf(\ T \L 1 /») , r=l-T/T c . (7) 

Hence, the normalized susceptibility x/L 2 ^ 11 becomes size-independent at t = and 
curves for different L plotted as functions of T exhibit a crossing point, similar to 
the behavior of the Binder cumulant. The collinear uud state of the triangular 
antiferromagnet breaks the Z3 lattice symmetry. The associated phase transition 
belongs to the universality class of the two-dimensional three-state Potts model, for 
which the critical exponents are known exactly including 77 = 4/15 125. . The bottom 
panel in figure [5] demonstrates that the normalized susceptibility curves for the above 
value of 77 cross precisely at the same point as the Binder cumulants. On one side this 
proves that the critical behavior for the transition into the uud state belongs to the 
three-state Potts universality class, on the other side it allows to employ stochastically 
less noisy Xq to determine the full transition boundary in the H—T plane. 

At low temperatures the transverse spin components exhibit a quasi-long-range 
order. The Binder cumulant method is not very convenient for a precise location of 
the corresponding Kosterlitz-Thouless (KT) transition. However, one can still use the 
normalized transverse susceptibility Xq/L 2 ^ 71 with the exact KT exponent 77 = 1/4 
[2"5] . We also measured the spin stiffness with respect to twists about the z-axis given 
for the triangular antiferromagnet by [22] 

P'= -^E«^ + ^SJ)) (8) 
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Here, dij — Vi — Vj and e stands for an arbitrary unit vector in the lattice plane. The 
spin stiffness shows a universal jump from zero to p s = 2Tkt/tt at the KT-transition. 
Measurements of p s (T) provide, therefore, an independent check for the susceptibility 
method. 

Results of various temperature and magnetic field scans are summarized in the 
H-T phase diagram presented in figure [3] The high-temperature paramagnetic phase 
remains invariant under the 0(2) x Z3 symmetry group, where the discrete symmetry 
Z3 is related to the permutation of three antiferromagnetic sublattices imposed by 
lattice translations: T a mQ — > e lC * a mQ. Three ordered phases in figure [3] break in 
different way the full symmetry group. The boundary between the paramagnetic and 
the uud state is formed by the second-order transition points, all other boundaries 
correspond to the KT-transitions. The notable feature of the phase diagram is a 
wide region with the collinear uud magnetic order [T5J [35] . At T = this collinear 
spin configuration appears only for one specific value of the applied magnetic field 
H = ^H s - Thermal fluctuations expand its stability range, such that the uud 
state has the highest transition temperature among the three ordered spin structures. 
The collinear ordering produces a clear 1/3-magnetization plateau, which is most 
pronounced in the plots dM/dH versus H. 

Overall the phase diagram shown in figure [3] is similar to the diagram proposed by 
Kawamura and Miyashita in |15j . The only qualitative difference concerns the low-field 
part. For temperature scans in a constant field H < \H S we always find two successive 
transitions: one from the paramagnetic state into the uud phase with the Z3-order 
of ^-components and, then, into the Y-state with an algebraic order of the transverse 
spin components. There is no direct transition between the paramagnetic state and 
the Y-phase suggested in [15] . The double transition is completely natural from the 
symmetry point of view because ordering processes for longitudinal and transverse spin 
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components are generally decoupled from each other. Note also that the Y-state and 
the mid-state should be separated from the disordered state in zero field by additional 
line(s) of phase transitions at H/J < 0.1. The present study remains uncertain about 
precise location of these transition boundaries. Numerical investigation of the low-field 
regime requires Monte Carlo simulations of substantially bigger clusters with L ~ 10 3 
because of the rapidly growing correlation length £ below T* ss 0.3 J [26] . 



3. Kagome antiferromagnet 

3.1. Multipolar phases in zero and finite magnetic fields 

The thermal order by disorder effect in the kagome antiferromagnet at zero field 
leads to a selection of coplanar states below the crossover temperature T* / J « 0.005 
[2TJI [21] . A true long-range order is not possible in 2D, though one can still discuss 
an asymptotic breaking of the spin-rotation symmetry at distances smaller than the 
correlation length £ ~ exp( J/T). The symmetry breaking in the ensemble of coplanar 
states is described by an unconventional octupolar order parameter [24] . For classical 
spin models magnetic octupoles are represented by the on-site third-rank tensor 

= S?S?S7 - l -Sf5 &1 - \s^5 ai - \s]5 a , . (9) 

The other common unconventional order parameter is the quadrupolar tensor [27] 

Qf = S?S? - l -5 aP . (10) 

Spin tensors ([9]) and ([T0|) transform under I = 3 and I = 2 representations of the 
0(3) rotation group, respectively. The usual (anti)ferromagnetic order parameter m q 
corresponds to the I = 1 representation. In a typical coplanar ground state of the 
kagome antiferromagnet spin on a given site may be oriented parallel to any of the 
three principal directions of the 120° spin triad. Thus, (Si) — and the lowest order 
spherical tensor, which captures the D3 (triatic) symmetry of the ensemble of coplanar 
states, is the uniform octupole tensor 

i 

The quadrupole tensor 

Q a0 = ^Y,^f) ( 12 ) 

i 

is the secondary order parameter and is induced via the rotationally invariant coupling 
term q^t^T^" . 

In an applied magnetic field the spin symmetry is reduced to the 0(2) group. Its 
irreducible representations are labeled by the projection of the angular momentum l z 
on the field direction. Therefore, different components of the same spin tensor may 
now describe different instabilities of the paramagnetic phase: 

Z, = ±l, Q za ,T zza , (13a) 

l z = ±2, Q ali ,T zaf> , (136) 

l z = ±3, T a ^ , a,/3, 7 = z,y. (13c) 

One can arbitrarily select either Q a @ or T^ 1 to characterize states with \l z \ = 1, 2. In 
the presence of a uniform magnetization m z the two tensors are coupled via a bilinear 




Figure 4. Magnetic field phase diagram of the Heisenberg kagome-lattice 
antifcrromagnet. Transition points determined by the Monte Carlo simulations 
are shown by squares. Solid lines are guides for the eye. 



term m z Q a ^T za ^ and become completely equivalent from the symmetry point of view, 
see also [25] , 

The two coplanar states stabilized in the kagome antiferromagnet at finite fields 
break differently the 0(2) rotation symmetry. The translationally disordered Y-state 
has the residual nematic symmetry C'2 and its order parameter is given by a 2 x 2 
symmetric traceless matrix (|I3bj) . The V-state may break the translational symmetry, 
then, its order parameter is the antiferromagnetic vector of the -\/3 x -\/3 structure |16j . 
If fluctuations are not able to select this periodic pattern, the proper order parameter 
of the V-state is transverse vector components of the quadrupolar tensor (|13a|) . An 
interesting observation is that the symmetry of the translationally disordered V-state 
coincides with the symmetry of the transverse ferromagnetic polarization m 1 allowing 
the interaction term m z Q za m a . For the Heisenberg model the classical constraint <j3j) 
ensures that = so that the corresponding coupling constant exactly vanishes. 
Nevertheless, the ground-state constraint may be modified, for example, by a single- 
ion anisotropy Dj^i^i' 2 - We have checked that in the easy-axis case D < the 
anisotropy produces a substantial transverse magnetization rrr 1 cx \D\ in the V-state 
of the kagome antiferromagnet. 

Partial breaking of the spin-rotation symmetry by multipolar order parameters 
leads, generally, to fractional topological defects [21]. The C2 symmetry of the Y-state 
is compatible with elementary half-quantum vortices, whereas the V-multipolar state 
supports only usual integer quantized defects. The universal jump of the spin stiffness 
at the KT-transition is given by [25, 29 



where v is the vortex winding number. For the half-quantum vortices with v — 1/2 
the jump in the spin-stiffness amounts to p s /T — 8/tt. 



3.2. Monte Carlo results 

A Monte Carlo investigation of the kagome antiferromagnet in external magnetic field 
has been performed using the same simulation protocol as described previously for 




Figure 5. Spin stiffness of the kagome antiferromagnet for two values of the 
applied magnetic field. Points are numerical data for the L = 36 cluster. Vertical 
arrows indicate positions of the Kosterlitz-Thouless transition obtained from the 
susceptibility data: T KT /J 6.5 X 10 -4 for H/J = 1 and T KT /J ss 9.3 X 10 -4 
for H/J = 5. Dashed lines represent the universal spin stiffness jump (f4) for 
v = 1/2 and v = 1 . 



the triangular antiferromagnet in section [2] To obtain a good statistics, one needs, 
however, to increase the number independent cooling runs to 30-50. We have studied 
periodic lattices with N = 3L 2 sites and L in the range 18-96. The obtained phase 
diagram of the kagome antiferromagnet is presented in figure 01 The Kosterlitz- 
Thouless transition boundaries for the two coplanar states were determined from the 
plots of the rescaled quadrupolar susceptibility with the exact KT-exponent 77 = 1/4 
[25] . We also measured the spin stiffness of the kagome antiferromagnet, which is 
obtained from the previous expression (JSJ by applying the renormalization factor 3 /4 
[24] . The temperature variation of p s for two values of the external field is shown in 
figure [5] The stiffness jump at the transition into the Y-state (H/J = 1) is equal 
with good accuracy to 8T/ir. This large jump demonstrates that the corresponding 
KT-transition is mediated by the unbinding of the half-quantum vortices. The jump 
of p s at H/J = 5 is substantially smaller and agrees with the standard value 2T/tt 
expected for the XY vector order parameter. 

Fractional vortices are incompatible with a conventional antifcrromagnetic 
ordering, therefore, the spin-stiffness data prove presence of the quasi-long-range 
nematic order in the low-field Y-state. In order to identify the type of symmetry 
breaking in the high- field V-state we consider in figure [5] the temperature evolution at 
H/J = 3 of the vector component of the quadrupole tensor Q 2 ± = (Q zx ) 2 + (Q zy ) 2 
and the transverse antifcrromagnetic amplitude toqj_ of the -\/3 x V3 configuration. 
The order parameters take on a finite value in the presence of a long-range order, but 
scale down as 1/iV, if correlations decay exponentially. The quadrupole tensor Q 2 ± 
demonstrates a very weak decrease with the lattice size at low temperatures, which 
is completely consistent with the power-law decay (QiQj) ~ and T-dependent 
r\ < 1/4. The antiferromagnetic order parameter goes down much faster, though still 
slower than 1/N. This is related to emergence of the Coulomb- type correlations that 
are typical to the constrained spin models, see, for example, [3D]- Overall, the Monte 
Carlo data show the absence of the translational symmetry breaking in the V-state. 

Thermal fluctuations have a strong effect in fields around ^H s , where they 
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Figure 6. Temperature dependence of the quadrupolc tensor , (upper panel) 
and the transverse antiferromagnetic order parameter of the v^3 X \/3 structure 
(lower panel) at H/J = 3. The vertical arrow indicates the Kosterlitz-Thoulcss 
transition temperature at Tkt / J ~ 6.7 X 10 — 3 . 

stabilize the collinear uud configuration between the two coplanar phases, see figured 
The corresponding 1/3-magnetization plateau survives to relatively high temperatures 
T/J ~ 0.05-0.2 [16]. The extreme stability of the collinear states is explained by the 
fact that they possess the largest possible number of zero-energy modes. Nevertheless, 
there is no symmetry breaking related to this selection and a collinear spin-liquid phase 
at the 1/3-plateau is connected to the paramagnetic phase by a simple crossover. 

Indeed, the collinear configurations preserve 0(2) rotations about the field 
direction and the only type of symmetry breaking at the plateau may be related to 
selection of a specific periodic pattern. The harmonic spectra of classical fluctuations 
are identical for different uud structures, though the anharmonic corrections are, 
generally, not |16) . The two particularly simple spin structures are the q = 
configuration and the -\/3 x state. In the dimer language, the two states correspond 
to the columnar array of dimers and the staggered arrangement, which maximizes the 
number of three-dimer hexagons. In our Monte Carlo simulations we did not find 
any significant tendency for the q — ordering in the classical model, though this 
configuration is favored by the zero-point fluctuations in the large-S* quantum kagome 
antiferromagnet [31 j - The numerical data for the antiferromagnetic order parameter 
of the longitudinal \/3 x y/3 spin structure are presented in figure [7J One finds again 
a fast, but slower than 1/N decrease of TOq with increasing lattice size. There is no 
evidence of the long-range ordering down to T/J ~5x 10~ 4 . 

The overall behaviour of the antiferromagnetic order parameter is consistent with 
the pseudo-dipolar correlations in the 2D Coulomb phase [3U] • At finite temperatures 
the Coulomb gas description of the ensemble of the collinear ground states at H = \H S 
is extended by allowing a nonzero density of monomers, which correspond to triangles 
with broken constraint condition The monomer doping destroys the power-law 
correlations at large distances producing a finite correlation length [32]. A detailed 
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Figure 7. Temperature variation of the longitudinal antiferromagnetic order 
parameter of the \/3 X V3 structure at the 1/3-magnetization plateau, H/J = 2. 

investigation of the crossover from short-distance power-law decay to exponential 
decrease at large distances requires systematic investigation of much bigger lattices. 

4. Summary and discussion 

In this paper, we have considered the effect of thermal fluctuations in external magnetic 
field for two classical frustrated antiferromagnets on 2D triangular and kagome lattices. 
The common property of the two models is selection of two coplanar and one collinear 
spin configurations by short-range fluctuations. For the triangular antiferromagnet 
such an order by disorder mechanism produces three ordered magnetic structures, see 
figure [3] The most dramatic effect concerns the collinear up-up-down configuration, 
which changes from being marginally stable at the mean-field level to having the 
highest transition temperature among the three states. 

For the kagome antiferromagnet, in addition to an arbitrary choice of the spin 
triad subject to the constraint ([3]) there is also a massive degeneracy related to different 
arrangements of the triad over the whole lattice. No periodic magnetic structure is 
selected down to at least T/J ~ 4 x 10~ 4 . We find instead unconventional magnetic 
ordering described by multipolar (tensor) order parameters, see the phase diagram of 
the kagome antiferromagnet in figure U] The most significant change again concerns 
the uud state, which does not break any symmetry and corresponds to a collinear spin 
liquid phase. Despite producing a pronounced 1/3-magnetization plateau this state 
is connected to a high-temperature paramagnetic phase by a crossover rather than a 
transition. 

The surprising feature of the phase diagram of the kagome antiferromagnet, 
figure 01 is an order of magnitude difference between transition temperatures of the 
high-field and the low-field multipolar state. Such a remarkable dissimilarity of the two 
temperature scales is produced by a combined effect of (i) 20% difference in the ground- 
state entropy, see ([2]) and ([4]), and (ii) different harmonic spectra with additional 
zero-energy modes present for the V-state [TB]. One can speculate, therefore, that if 
other terms are added to the spin Hamiltonian (JTJ , the system would more easily find 
an optimal periodic structure in the V-state rather than in the more degenerate Y- 
state. Let us also note, that there is an interesting similarity between the phase 
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diagram in figure U and the experimental phase diagram of the 3D hyperkagome 
antiferromagnet Gd3Ga 5 0i2, which orders only in finite magnetic fields [33]. Thus, 
a Monte Carlo investigation of the magnetic diagram of the hyperkagome model 
is a pressing extension of the present study. One should be also intersted in the 
effect of the single-ion anisotropy on the phase diagram of the triangular Heisenberg 
antiferromagent in relation to the experimental studies of RbFe(Mo04)2 [13 HB] and 
Rb 4 Mn(Mo0 4 )3 [H]- 
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